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Abstract. Modelling global disc galaxies is a difficult task which has previ- 
ously resulted in the small scale physics of the interstellar medium being greatly 
simplified. In this talk, I compare simulations of galaxies with different ISM 

^ \ properties to determine the importance of the ISM structure in the star forma- 

O I tion properties of the disc. 



Introduction 



JL ' One of the main problems in understanding galaxy formation is that star forma- 

5^ . tion is extremely complex. The gas out of which stars form is a turbulent mix 

j/3 I of forces in which gravitational collapse, thermal pressure, magnetic fi elds, cos- 



C^ . mic ray pressure and energy from supernovae all fight for dominance (JMacLow 

120041 1. 

This leaves galaxy simulators with a choice; either to model the interstellar 



^ \ gas in detail but restrict their study to a small patch of the galaxy (e.g. ISlvz et al.. 

^ ■ I2OO5I ) or t o simulate the entire galaxy but use a greatly simplified model for the 



ISM (e.g. iLi et al.ll2005l ). The former approach allows the inclusion of many 



m 

(^ ■ more of the important physical processes but is unable to tell us anything about 

the global evolution of the disc. The latter, meanwhile, reveals properties of the 
whole galaxy including star formation histories and global structures, but it is 

C^^ \ impossible to judge the effect the simple ISM model is having on these results. 

^^ ' However, recent numerical simulations are now able to include a more 

complex multiphase IS M in these global models ( Tasker &: Brvarj |2007| . l2006l : 

'%_] . in a more realistic way in galaxy formation. This has particular baring on cos- 

C^ \ mological simulations, where the resolution of the ISM of individual galaxies is 

still beyond our reach. 

2. Numerical method 



Using the AMR code, Enzo ( Brvan &: Normanlll997l ). we compared three models 



of isolated galaxy discs where we varied the properties of the interstellar gas. 
In all cases, the simulations started with a rotating Milky Way-sized exponen- 
tial disc of gas sitting in a static NFW dark matter potential. The set-up is 
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Figure 1. Gas density projections at 388Myrs for (left to right) ISM #1, 
ISM #2 and ISM #3, all without feedback. Scale is to the base-10 logarithm 
with units M0Mpc~^ and each image is 60kpc across. 



described in more detail in iTasker &: BrvanI ( 20061 . [2003) • In our first galaxy 
model (ISM #1), the gas was allowed to radiatively cool to 300 K. In our second 
model (ISM #2), an additional background photoelectric heating term was in- 
cluded while in our last disc (ISM #3), the ISM was kept at a fixed temperature 
of 10,000 K, in keeping with previous simulations. All simulations included star 
formation and ISM models 1 and 2 also included feedback from Type II SNe 
(this was impossible to include for the third, isothermal, disc). 



3. The structure of the ISM 



Figure [T] shows projections of the gas density of each galaxy disc 377Myrs after 
the start of the simulation. At this time, the disc has undergone its initial 
fragmentation (and subsequent star burst) and the evolution is now slow. In 
Figure [H all images are for runs without stellar feedback. With our first ISM 
model, shown on the far left, we see that the gas has fragmented out to a well- 
defined radius. Beyond this point, there is still a significant amount of gas, 
but it is stable to gravitational fragmentation and does not collapse to form 
stars. When we include photoelectric heating in ISM #2 (middle image) we see 
a notably different gas structure. In particular, there exists large voids of hot, 
low density gas. The porous nature of this ISM has been observed both in our 
own galaxy and (perhaps most dramati cally) in the HI ma p of the LMC. This 
result agrees with recent simulations of IWada et al.l ( 2000! ). who suggest these 
holes are not the results of SNe remnants, but rather the product of gravitational 
and thermal instabilities. Our final ISM model, where the temperature is fixed, 
shows another distinct structure. Due to the fixed high temperature of 10^ K, 
the Jean's Length is higher than in the other discs, causing the star-forming 
knots of gas to be much higher in mass. This has the effect of producing very 
large star clusters that are confined to the densest, inner regions of the disc. 
Their formation also produces voids in the gas distribution, but this is due to 
gas deficit, not to instabilities. 

The pressure distribution of these discs is also interesting. Without feed- 
back, discs with ISM ^1 and ISM #2 are largely in pressure equilibrium, as 
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Figure 2. PDF of the volume weighted gas density at 377 Myrs for the three 
ISM types (running left to right). A log normal fit is shown as a solid line in 
the first panel. 



predicted by the analytical models of lMcKee &: Qstriken ( 19771 ). The isothermal 
disc cannot be, since fixing the temperature requires the pressure to be propor- 
tional to the gas. When feedback is introduced, this situation changes. SNe 
energy causing streams of gas to be blown both in the disc's plane and off its 
surface, creating a galactic fountain. 

The structure of these ISMs can be examined more quantitatively by looking 
at their ID PDFs, as shown in Figure [2j Comparing all three of the ISMs across 
the panels, we see a significant amount of substructure in the low density gas, but 
the high density end of the PDFs are surprisingly uniform. At these densities, 
all profiles are well fitted by a lognormal distribution. This remains true even 
when feedback is included. 



4. Observational comparison 

The left-hand plots in Figure [3] shows the star formation history for each galaxy 
disc. Without feedback, our ISM #1 disc converts all the available gas into 
stars, halting further star formation. By contrast, the addition of background 
heating quenches star formation by raising the temperature of the coldest gas 
and allowing the disc to show the beginning of self-regulation. The addition of 
feedback, however, is a much stronger effect, with the added energy destroying 
the star forming knots and thereby increasing the available gas at later times. 
The isothermal disc also shows a flattening in the star formation rate, but this is 
likely due to the confinement of the star formation to the denser parts of the disc 
which slows down the rate of gas consumption. We can also cornpare our results 
with the widely observed Kennicutt-Schmidt law ( KennicuttI Il989l ). as shown 
in the right-hand plot of Figure El Here, we see that the observed gradient is 
reproduced well in both the first two ISM types, but less well in our isothermal 
disc. We do however, persistently overestimate the rate of star formation, even 
when feedback is included. This is likely due to our star formation recipe and 
is something to address in later work. 
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Figure 3. The star formation history (left) and global Schmidt-Kcnnicutt 
law (right) for ISM models 1, 2 and 3 (top to bottom). 



Conclusions 



One of the most surprising results in this research is that despite significant 
structural differences in the ISM, the star formation properties remain largely 
unaffected. The reason for this can be seen in Figure [21 where we see that the 
biggest differences between our three ISM types occur in the low to medium 
density gas. The high density, star-forming gas, meanwhile, forms a consistent 
lognormal profile in all cases. This result is largely positive; it suggests that 
a decent subgrid model for star formation can be used where a detailed ISM 
model is not possible, such as large-scale cosmological simulations. However, 
it also implies that a through understanding of star-formation cannot be found 
from the Kennicutt-Schmidt law, which is l argely insensitive to the input physics. 
This work is presented in greater detail in lTasker &: BrvanI ( 20071 ). 
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